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We revisit the electroweak phase transition and the critical bubble in the scale invariant two 
Higgs doublet model in the light of recent LHC data. Moreover, the sphaleron decoupling condition 
is newly evaluated in this model. The analysis is done by using the resumed finite-temperature 
one-loop effective potential. It is found that the 125 GeV Higgs boson inevitably leads to the strong 
first-order electroweak phase transition, and the strength of which is always large enough to satisfy 
the sphaleron decoupling condition, vn/Tn > 1.2, where Tjv denotes a nucleation temperature and 
wjv is the Higgs vacuum expectation value at Tjv. In this model, even if the Higgs boson couplings 
to gauge bosons and fermions are similar to the standard model values, the signal strength of the 
Higgs decay to two photons is reduced by 10% and the triple Higgs boson coupling is enhanced by 
82% compared to the standard model prediction. 


I. INTRODUCTION 

One of the observational facts that needs new physics 
beyond the standard model (^M) is the baryon asymme¬ 
try of the Universe (BAU) [y, 

— = (8.59 ±0.11) X 10-^1 (Planck) (1) 
s 

where ns (s) denotes the baryon number (entropy) den¬ 
sity. Although many mechanisms that can explain the 
observed value exist in the literature, electroweak baryo- 
genesis (EWBG) Q is the only scenario that is ripe for 
verification by collider experiments, such as the Large 
Hadron Collider (LHC), and by low energy experiments, 
such as the electric dipole moments of the neutron, atoms 
and molecules. Since EWBC is intimately connected 
to Higgs physics, the establishment of the Higgs sec¬ 
tor plays an essential role in testing it, and the dis¬ 
covery of the Higgs boson at the LHC in 2012 Q 
is the first step toward the collider probe of EWBC. 
Indeed, since the Higgs boson mass that is one of the 
relevant parameters for the electroweak phase transi¬ 
tion (EWPT) has been measured with 0.2% accuracy, 
rriH = 125.09 ± 0.21 (stat.) ± 0.11 (syst.) CeV the 
feasible regions of EWBC have been narrowed down in 
various models @ . In upcoming experiments, such as 
the LHC Run-H and the High-Luminosity LHC Q, the 
Higgs boson couplings to the SM particles would be mea¬ 
sured with better precision, and the international linear 
collider (ILC) Q has the great capability of measuring 
the triple Higgs boson coupling, which may yield a deci¬ 
sive clue to the EWBC hypothesis. 

In order for EWBC to be successful, the EWPT has 
to be strongly first order. The properties of the EWPT 
are related not only to the Higgs boson mass and model 
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parameters but also to electroweak symmetry breaking 
mechanisms. An interesting possibility is the so-called 
Coleman-Weinberg (CW) mechanism I, 0 in which 
quantum effects induce the symmetry breaking. The 
scale invariant two Higgs doublet model (SI-2HDM) [n^ 

a is one of such examples, 0 and the previous work 
shows that the SI-2HDM can have the strong first- 
order EWPT. At the time of their analysis, however, 
the masses of the Higgs boson and top quark were not 
known. Moreover, on the theoretical front, neither a ther¬ 
mal resummation for the effective potential nor the eval¬ 
uation of a baryon number preservation condition (also 
called a sphaleron decoupling condition) were conducted 
in Ref. 

In this Letter, we update the analysis of the EWPT in¬ 
cluding the evaluation of bubble wall profiles, and obtain 
the sphaleron decoupling condition by taking the recent 
LHC data into account. In our study, we use the finite- 
temperature one-loop effective potential with daisy re¬ 
summation. The phenomenological consequences of the 
sphaleron decoupling condition is also briefly discussed. 
As studied in the previous works [ 13,13 , we evaluate the 
deviations of the Higgs boson couplings from their SM 
values in the region where the strong first-order EWPT 
is achieved. 

The paper is organized as follows. We give a quick re¬ 
view of the SI-2HDM in section m and the Higgs boson 
couplings are presented in section IHII The sphaleron de¬ 
coupling condition and the critical bubbles are discussed 
in section ITVl We show our results in section |Vl and con¬ 
clusions and discussions are given in section IVTl 


^ The CW mechanism does not work in the SM since the top quark 
is too massive to give the stable vacuum. 
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II. THE MODEL 


The ST2HDM is a minimal scale invariant extension 
of the SM by adding another Higgs doublet field. The 
most general Higgs potential at the renormalizable level 
is given by 

+ A4($1$2)($J$i) + I 

+ A 6 ($J$i)(d>l$ 2 ) + A 7 ($J$ 2 )(«'I$ 2 ) + h.C.j. 

( 2 ) 


After two Higgs doublets get vacuum expectation values 
(VEVs), they are cast into the form 


$i(x) = 




<Pt (a^) 

- hi{x) - 


^{x)) 


^ = l,2, (3) 


where vi = vcosP and V 2 = usin/? with 0 < /3 < tt/2, 
and V ~ 246 GeV. In order to avoid Higgs-mediated flavor 
changing neutral current (FCNC) processes at the tree 
level, we impose a Z 2 symmetry (d’l —^ —$ 1 , $2 —>■ $ 2 ), 
which leads to Ae = A 7 = 0 [1^. The phase of A 5 is 
removed by an appropriate field redefinition of the Higgs 
doublets, so that CP is conserved. 

Following a method by Gildener and Weinberg [13 , we 
consider the FW symmetry breaking in a flat direction. 
The tree-level effective potential takes the form 

T/ / \ "^1 4 I '^2 4 A345 22 / A \ 

Ho(</ 51,(/32) = Y'f'l + ■§■‘^2 + —(4) 


where ifi and ip 2 are the constant background fields of 
the two Higgs doublets. 

The tadpole conditions that are defined as the first 
derivatives of Vq with respect to (pi ^2 give the following 
conditions: 


A345 + \/AiA2 — 0 , Aii:3 — X2V2, 


( 5 ) 


where A 345 = A 3 -I- A 4 -I- A 5 . With these conditions, it 
follows that Vo(vi,V 2 ) = 0. Moreover, since the mass 
matrix of hi and /12 is written as 


MLe 


f Xivl A345?^l'C2 
\A345r!i'i;2 A2r’| 


( 6 ) 


one finds det(A4tree) = 0 using Eq. ([5]). The appearance 
of the massless particle is the consequence of the classical 
scale invariance. We define h and H as the mass eigen¬ 
states of the CP-even Higgs bosons, which are obtained 
by 


where —-Kj^ < a < 0. In the following discussion, h 
is the SM-like Higgs boson whose mass is zero at the 
tree level and is generated by the quantum corrections. 
It can be proved that a = P — nf 2 at the tree level, 
and consequently, the Higgs boson couplings to the gauge 
bosons and fermions are the same as those in the SM. 

As mentioned above, h becomes massive as the result 
of the radiative EW wmmetry breaking. The one-loop 
effective potential is [ 3 , 1 ^ 



where ip\ and i = H, A, H^, W^, Z, t, b, 

and Ci = 3/2 (5/6) for scalars and fermions (gauge 
bosons) and /2 denotes a renormalization scale. A and 
are the physical CP-odd and charged Higgs bosons, 
respectively, rii are the degrees of freedom and the statis¬ 
tics of the particle i: 

77/7 — 77j4 — 1, 77.774 — 2 , 774^4 — 3 ■ 2, 

nz = 3, nt = Ub = -12. (9) 

The field-dependent masses can be written as 777 ^ = 
where rrii are the corresponding masses in the 
vacuum, so that Vi {(p) is reduced to 


Vi{ip) = Atp"^ -I- Bif'^ log 


‘fl 

-9 ’ 
/i 2 


( 10 ) 


with 


A = 


m. 





( 11 ) 


As can be seen from the tadpole condition of Vi ((/?), we 
have a relationship between the scale of v and the renor¬ 
malization scale / 2 , i.e., as the con¬ 

sequence of dimensional transmutation. From Eqs. m 
and (HU, it is easily checked that the vacuum energy 
becomes Vi( 7 ;) = —Bv^/2, which implies that the elec- 
troweak symmetry is broken unless B is negative. It 
should be noted that since A and B are the same order in 
the coupling, i.e., 0{g'^), where g collectively denotes the 
coupling constants in this model, A/B should be 0(1), 
so the \og{v‘^/p?) ~ 0(1). In other directions, however, 
A may be 0{g^) and thus log{v^/fl^) ~ l/g^, which may 
invalidate the perturbative calculation, as advocated in 
Ref. [13]. 

The mass of h is obtained by taking the second deriva¬ 
tive of Viipp) and evaluating it at (/? = ri. 


ml = 


d'^Viip) 

dip'^ 


^p—v 


8Bv^. 


( 12 ) 



/COS a — sin a 
\^sina cos a 



We remark that thanks to the loop contribntions from 
i7, A and B can be positive in contrast to the SM 
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case, rendering positive. Interestingly, once rrih = 
125 GeV is fixed, the possible ranges of rriH, mA and 
mH± are restricted. In this Letter, we consider a case 
that niA = rriHi in order to satisfy the constraint coming 
from the p parameter [2lj. Therefore, the heavy Higgs 
mass scales are specified by only two parameters. In what 
follows, ruH and niA are chosen. 


III. HIGGS BOSON COUPLINGS 


The Higgs boson couplings to gauge bosons and 
fermions normalized to the SM values are, respectively, 
given by 


Kv 


„SI-2HDM 

yhvv 


9hvv 




„SI-2HDM 

9hff 


„SM 

9hff 


(13) 


where f = u,d, 1. As discussed in section H, Ky = Kf = 1 
due to a = /3 — it 12 at the tree level. Even in such a sit¬ 
uation, the so-called nondecoupling effects may appear 
in the loop processes. For instance, as pointed out in 
Ref. [ 2 ^, the h —> 77 mode may be significantly modi¬ 
fied by the charged Higgs boson loop. The Higgs signal 
strength of h —>• 77 is defined as 


M 77 — 


a{pp —>• h)si-2HDMBr(/i —>• 77 )si-2HDM 


1 -b 


a{pp /i)sMBr(/i 77 )sm 
Ah± ' 


Asm 


(14) 


where Asm = —6.49 and Ah± = —th± (l — 
'’■ff±/('''ff±)) with T^± = 4m^±/m^, and / is a loop 
function defined in Ref. (^ . 

The another nondecoupling effect may appear in the 
triple Higgs boson coupling. The deviation of the triple 
Higgs boson coupling from its SM value is defined as 

\ SI-2HDM _ \ SM 

A \ ^hhh ^hhh /ic:\ 

^-^hhh = -Tsm-■ 

'^hhh 

In this analysis, we use the following expression as the 
SM prediction [ 2 ^ 


i,SM _ 
^hhh — 


1 -b 


9ml 


E 

i=W,Z,t,b 


n^ 




(16) 


Note that the dominant one-loop contribution comes 
from the top quark loop, which renders Xuhh smaller com¬ 
pared to the leading result. In the ST2HDM, the triple 
Higgs boson coupling to leading order is simply expressed 
in terms of mh and v [ 2 ^ 


\SI-2HDM 

^hhh 


d^Vi{p) 

dif^ 


ip—V 


= 40Bv = 



V 


Unlike the ordinary 2HDM, the leading result in the ST 
2HDM does not same as the leading one in the SM even 
in the case that (3 — a = 7 r/ 2 , which reflects the different 
origins of the electroweak symmetry breaking. 


IV. SPHALERON DECOUPLING CONDITION 
AND CRITICAL BUBBLES 

In EWBG, in order to preserve the generated BAU un¬ 
til today, the sphaleron process must be decoupled right 
after the electroweak symmetry breaking. This condition 
(the so-called sphaleron decoupling condition) is given by 

rg^(r) <iJ(T), (18) 


where (T) is the baryon number changing rate in the 
broken phase, and H{T) is the Hubble constant. Eq. 
can be translated into 


v{T) 

T 


4tt£{T) 

= Csph('H), 


42.97-b logAT 


(19) 


where the sphaleron energy is denoted as Esph = 
4TTv(T)£(T)/g2, with 52 being the SU(2) gauge coupling. 
M represents the translational and rotational zero-mode 
factors of the fluctuations about the sphaleron. 

In our numerical analysis, we first evaluate both Tq 
and VC, where Tc stands for a critical temperature at 
which the two degenerate minima appear in the effective 
potential, and vc is the VEV of the Higgs fields at Tc- 
The EWPT is studied in the direction of p, and tan/3 
is fixed by that at T = 0. We use the resummed finite- 
temperature one-loop effective potential 


Ves{^,T) 


647r2 









( 20 ) 


where 

poo 

lB,F{a^) = J dx log [1 =F j , ( 21 ) 

with the upper (lower) sign for bosons (fermions). 

are the thermally corrected boson masses de¬ 
fined as Mf{ip,T) = + ni(T) where ni(r) are 

the finite-temperature self-energy. Here, we consider the 
leading 0{T^) terms 


n$(T) = 


T2 

1^2 [ 


6 m^ 


3m7 




m\ 


2w? 


H± 


+ -b ml) 


Hw (T) =2gir^ , Hb (T) = 2glT^ 


( 22 ) 

(23) 


where H^ for the Higgs bosons, Hw and Hb for the SU(2) 
and U(l) gauge bosons, respectively. Note that the only 
longitudinal part of the gauge boson self-energy is ther¬ 
mally corrected. 

After finding Tc, we evaluate the sphaleron energy at 
that temperature (for a detailed calculation, see, e.g., 
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Refs. Since the dominant contribution in the 

right-handed side of Eq. m comes from £(T), we ne¬ 
glect the logarithmic terms in our numerical analysis. 

It should be noted that the EWPT does not start at 
T = Tc- In order for the EWPT to occur, the bubbles 
must be nucleated at somewhat below Tc- Only bubble 
that has some critical size, which is called the critical 
bubble, can grow. The EWPT proceeds to develop if 
the bubble nucleation rate is larger than a certain value, 
and then the Universe is finally filled with the broken 
phase. We define the nucleation temperature (Tjv) by 
the condition 


TN{TN)/H^iTN) = HiTN), (24) 


where T]s[(Tn) denotes the bubble nucleation rate per 
unit time per unit volume at Tn [HI. Prom Eq. (|24ll . it 
follows that 

Ecb(TAf) _ 3 £’cb(TAf) 

rji rt O rji 

In ^ 

= 152.59 - 21ogg.(riv) - 4 log , (25) 

where Ecb(T/v) is the energy of the critical bubble and 
S'* (Tat) represents the degrees of freedom of the relativis¬ 
tic particles at T^. As seen from Eq. (ESI) . Ech/T < 150 
is needed for development of the EWPT. 

We closely follow a method in [s^ to evaluate Ech{T). 
The critical bubbles are estimated from the following en¬ 
ergy functional 


Ecb(T) 






FIG. 1. Contours of the Higgs boson mass and vcjTc in 
the {niH, rnA) plane. The solid line in black corresponds to 
nih = 125 GeV, and the dashed lines in gray indicate nih = 
80, 200 and 250 GeV from bottom to top. The each contour 
in white represents vc jTc = 1.1, 1.5, 2.0, 3.0 and 4.0 from 
top to bottom. 


with the boundary conditions: dhi^2{0/d^\^=o ~ ^ 
/ii, 2 (C = oo) = 0. With those solutions, Ech{T) is evalu¬ 
ated. 

It is known that the bubble solutions are approxi¬ 
mately given by a kink-configuration 


+ Ves{<^i,^2,T) , 


(26) 


Pi{r) ^ ViiT) 


1 — tanh 


r — R 


(30) 


where the gauge fields are assumed to take the pure-gauge 
configuration so that they do not contribute to the en¬ 
ergy of the critical bubbles. The classical Higgs fields are 
parametrized as 

^ 72 Ciw) ’ ^ 72 C2(r)) ’ 

where pi{r) = p{r) cos/3, P 2 (j') = pir) sin (3, r = \x\, and 
here tan /3 is fixed by that at T = 0 as mentioned above. 
In the numerical analysis, it is convenient to change r 
and Pi into the following dimensionless quantities: 


where R and are the radius and wall width of the 
bubbles, respectively. We use this as the initial config¬ 
uration to derive the bubble solutions by using the re¬ 
laxation method. Eor more details about the numerical 
method, see, e.g.. Ref. [30l| . 

V. NUMERICAL RESULTS 

In the SI-2HDM, there are five parameters in the tree- 
level potential: 


e = u(T)r, hiiO 


Pijr) 

v{T) cos/3’ 


h2{0 


P2{r) 


v{T) sin/3 
(28) 


Ai, A 2 , A 3 , A 4 , A 5 . (31) 

In our analysis, we replace them with the following phys¬ 
ical parameters: 


The profiles of hi (^) are obtained by solving the equations 
of motion 

_ 72 ^^1(2) \__ dVeS _ 

7 dC V d^ J u"‘(T) cos^/3(sin^/3) dhi( 2 ) ’ 

(29) 


niH, niA, rnH±, /3, v. (32) 

We take mH± = rriA as mentioned in sectionlHl Since VeS 
does not depend on tan P explicitly, the results obtaining 
from it do not either, except for the cutoff of the model, 
as will be discussed in the following. 
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FIG. 2. Bubble profiles of /ii, 2 (C) at T = Tn- In this plot, we 
set TTih = 125 GeV and mn ~ rnA = 382 GeV. 


In Fig. [1] we show the contours of the Higgs boson 
mass and vciTc in the (to//, itia) plane. The black 
solid (gray dashed) line indicates the parameter region 
where to/j = 125 (80, 200, 250) GeV. As can be seen 
from Eq. ca), the Higgs boson mass gets larger as to// 
and niA increase. The white contours represent the mag¬ 
nitude of VC /Tc- These contours indicate that the size of 
vc/Tc becomes smaller as to // and uia get heavier. Since 
the thermal effects from the heavy Higgs bosons cause the 
first-order EWPT in this model, vc jTc would be propor¬ 
tional to V ^ ^. mf, i = H,A, from the high- 

temperature approximation argument. This may explain 
the behavior of vcjTc qualitatively. As a benchmark 
point, we take rtih = 125 GeV and to// = toa = 382 GeV. 
In this case, we find that vcITq = 211 GeV/91.5 GeV = 
2.31 and Csph(T’c) = 1.23. Therefore, even though ^gph is 
greater than the conventional criterion by 23%, vcjTc 
is large enough to satisfy the sphaleron decoupling con¬ 
dition. It is also found that Csph(Tc') is almost constant 
on the black line while vcjTc gets slightly weaken in the 
region where to// > 500 GeV since the thermal effect of 
H loop is suppressed. 

In Fig. [21 the profiles of hi^2{0 is shown, here hi{^) = 
by construction. Our numerical calculation shows 
that vn/Tn = 229 GeV/77.8 GeV = 2.94, Csph(TAr) = 
1.20 and Acb(7Ar)/T/v = 151.7. The degrees of the super¬ 
cooling is about 15%, i.e., {Tq — Ta[)/Tc = 0.15, which 
is more or less the same as the previous estimate [T^ . 

Let us briefly make a comparison between the SI- 
2HDM and the minimal supersymmetric SM (MSSM). 
Unlike the SI-2HDM, the supercooling in the MSSM case 
is rather small, e.g., 0(10“^) [1^, and the bubble wall 
width in the SI-2HDM is thinner than that in the MSSM, 
which is due to the stronger first-order EWPT compared 


ruH 

382 GeV 

vc/Tc 

211 GeV/91.5 GeV = 2.31 

Csph(T’c) 

1.23 

vn /Tn 

229 GeV/77.8 GeV = 2.94 

Csph(7jv) 

1.20 

Ea,(Tn)/Tm 

151.7 

Ky 

1.0 

Kf 

1.0 

M 77 

0.90 

AA/i/i/i 

82.1% 

A 

6.3 TeV 


TABLE I. The benchmark point for the strong first-order 
EWPT and the nucleation of the bubbles, where we take 
nih = 125 GeV and mn = toa = m[f±. For the evaluation of 
A, tan /3 = 1 is used. 


to the MSSM case. Since a CP violating source term may 
be proportional to the gradient of the bubble wall , the 
baryon number generation may be more efficient than the 
MSSM case. To this end, of course, the current model has 
to be extended to have an extra source of CP violation. 

In Table HI our numerical results in a benchmark point 
are summarized. In the SI-2HDM, the strong first-order 
EWPT is the inevitable consequence from the require¬ 
ment of the 125 GeV Higgs boson. In this case, the sig¬ 
nificant deviations may appear in and Xhhh- We leave 
the study on the detectability of the heavy Higgs bosons 
to future work. 

Finally, we comment on the cutoff scale (A) of this 
model. Here, A is determined by a scale at which 
|Ai| > Itt is obtained. In doing so, we use the one-loop 
renormalization group equations [32|. As an example, 
tan /I = 1 is taken. It is found that A = 6.3 TeV, which 
is extremely low compared to a typical grand unihcation 
scale, ^ 10^® GeV. Since Ai oc tan^ (3 and A 2 oc 1/ tan^ /3, 
the cases for tan /3 > 1 and tan /3 < I yield the lower 
cutoff scales than 6.3 TeV. Our analysis has reconfirmed 
the previous results Bm. 


VI. CONCLUSIONS AND DISCUSSIONS 

We have revisited the EWPT and the profiles of the 
critical bubbles in the SI-2HDM in the light of the 125 
GeV Higgs boson. We improved these analyses by using 
the finite temperature one-loop effective potential with 
thermal resummation. In this model, the heavy Higgs 
mass scales are fixed to be consistent with the mh = 125 
GeV. In our benchmark point, to// = toa = to//± = 382 
GeV, we found that vc/Tc = 211 GeV/91.5 GeV = 2.31 
and Csph(Tc') = 1.23. At the nucleation temperature, 
they are changed into vn/T^i = 229 GeV/77.8 GeV = 
2.94 and Csph(7Ar) = 1.20. Even though ^gph in the SI- 
2HDM is greater than the conventional criterion by about 
20%, the first-order EWPT is strong enough to satisfy the 
sphaleron decoupling condition. 

We also studied the deviations of the Higgs boson cou- 
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plings from the SM predictions. It was found that even 
though the Higgs boson couplings to the gauge bosons 
and fermions are SM like, the significant deviations may 
appear in the —>■ 77 mode and the triple Higgs boson 
coupling due to the nondecoupling effects of the heavy 
Higgs boson loops. In our benchmark point, the Higgs 
signal strength of —>■ 77 is reduced by 10 % and the 
triple Higgs boson coupling is enhanced by 82.1%. Such 
deviations may be detectable in the future experiments 
such as the High-Luminosity LHC 0 and the ILC Q. 

There are some issues to be solved. In order to obtain 
the baryon asymmetry, an extra source of CP violation 
is needed as mentioned in the previous section. To this 
end, we may augment this model by adding the extra 
fermions in a scale-invariant way. The current analysis 
would not be much modified as long as the strength of the 


interactions between new particles and the Higgs boson 
are moderate. Furthermore, since the cutoff of the model 
is rather small, the UV completion is needed. However, 
constructing a complete model is beyond the scope of this 
Letter. 


ACKNOWLEDGMENTS 

The work of K.F. is supported by Research Fellow¬ 
ships of the Japan Society for the Promotion of Science 
for Young Scientists. E.S. is supported in part by the 
Ministry of Science and Technology of R. O. C. under 
Grant No. MOST 104-281 l-M-008-011. 


[1] P. A. R. Ade et al. [Planck Collaboration], Astron. As- 
trophys. 571, A16 (2014). 

[2] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposh- 
nikov, Phys. Lett. B 155 (1985) 36. For reviews on 
electroweak baryogenesis, see A. G. Cohen, D. B. Ka¬ 
plan and A. E. Nelson, Ann. Rev. Nucl. Part. Sci. 43 
(1993) 27; M. Quiros, Helv. Phys. Acta 67 (1994) 451; 
V. A. Rubakov and M. E. Shaposhnikov, Usp. Fiz. Nauk 
166 (1996) 493; K. Funakubo, Prog. Theor. Phys. 96 
(1996) 475; M. Trodden, Rev. Mod. Phys. 71 (1999) 
1463; W. Bernreuther, Lect. Notes Phys. 591 (2002) 
237; J. M. Cline, arXiv:hep-ph/0609145 ; D. E. Morris¬ 
sey and M. J. Ramsey-Musolf, New J. Phys. 14, 125003 
(2012); T. Konstandin, Phys. Usp. 56 (2013) 747 [Usp. 
Fiz. Nauk 183 (2013) 785]. 

[3] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 

1 ( 2012 ). 

[4] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 
716, 30 (2012). 

[5] G. Aad et al. [ATLAS and CMS Collaborations], 
arXiv:1503.07589 [hep-ex]. 

[6] T. Cohen, D. E. Morrissey and A. Pierce, Phys. Rev. D 
86, 013009 (2012); D. Curtin, P. Jaiswal and P. Meade, 
JHEP 1208, 005 (2012); M. Carena, G. Nardini, 
M. Quiros and G. E. M. Wagner, JHEP 1302, 001 
(2013); K. Krizka, A. Kumar and D. E. Morrissey, 
arXiv:1212.4856 [hep-ph]. 

[7] [GMS Gollaboration], arXiv: 1307.7135; [ATLAS Collab¬ 
oration], arXiv: 1307.7292 [hep-ex]. 

[8] H. Baer, T. Barklow, K. Fujii, Y. Gao, A. Hoang, 
S. Kanemura, J. List and H. E. Logan et al., 
arXiv:1306.6352 [hep-ph]. 

[9] S. R. Coleman and E. J. Weinberg, Phys. Rev. D7 (1973) 
1888. 

[10] E. Gildener and S. Weinberg, Phys. Rev. D13 (1976) 
3333. 

[11] K. Inoue, A. Kakuto and Y. Nakano, Prog. Theor. Phys. 
63 (1980) 234. 

[12] K. Funakubo, A. Kakuto and K. Takenaga, Prog. Theor. 
Phys. 91, 341 (1994). 

[13] K. Takenaga, Prog. Theor. Phys. 92 (1994) 987. 

[14] K. Takenaga, Prog. Theor. Phys. 95 (1996) 609. 


[15] J. S. Lee and A. Pilaftsis, Phys. Rev. D 86 (2012) 035004. 

[16] G. T. Hill, Phys. Rev. D 89 (2014) 7, 073003. 

[17] S. Kanemura, Y. Okada and E. Senaha, Phys. Lett. B 
606 (2005) 361. 

[18] K. Fuyuto and E. Senaha, Phys. Rev. D 90, no. 1, 015015 
(2014). 

[19] S. L. Glashow and S. Weinberg, Phys. Rev. D 15 (1977) 
1958. 

[20] R. Jackiw, Phys. Rev. D 9, 1686 (1974). 

[21] H. E. Haber and A. Pomarol, Phys. Lett. B 302 , 435 

(1993); A. Pomarol and R. Vega, Nucl. Phys. B 413, 
3 (1994); P. H. Ghankowski, M. Krawczyk and J. Zo- 
chowski, Eur. Phys. J. C 11, 661 (1999); J.-M. Ger¬ 

ard and M. Herquet, Phys. Rev. Lett. 98, 251802 
(2007); W. Grimus, L. Lavoura, O. M. Ogreid and 
P. Osland, Nucl. Phys. B 801, 81 (2008); S. de Visscher, 
J. M. Gerard, M. Herquet, V. Lemaitre and F. Maltoni, 
JHEP 0908, 042 (2009). 

[22] 1. F. Ginzburg, M. Krawczyk and P. Osland, In *Seog- 
wipo 2002, Linear colliders* 90-94 hep-ph/0211371 . 

[23] D. McKeen, M. Pospelov and A. Ritz, Phys. Rev. D 86 
(2012) 113004. 

[24] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, 
Front. Phys. 80 (2000) 1. 

[25] S. Kanemura, S. Kiyoura, Y. Okada, E. Senaha and 
G. P. Yuan, Phys. Lett. B 558 (2003) 157; S. Kane¬ 
mura, Y. Okada, E. Senaha and C.-P. Yuan, Phys. Rev. 
D 70 (2004) 115002. 

[26] D. Chway, T. H. Jung, H. D. Kim and R. Dermisek, Phys. 
Rev. Lett. 113 (2014) 5, 051801. 

[27] M. E. Carrington, Phys. Rev. D 45, 2933 (1992). 

[28] N. S. Manton, Phys. Rev. D 28, 2019 (1983). 

[29] F. R. Klinkhamer and N. S. Manton, Phys. Rev. D 30, 
2212 (1984). 

[30] K. Funakubo and E. Senaha, Phys. Rev. D 79 (2009) 
115024. 

[31] A. D. Linde, Nucl. Phys. B 216, 421 (1983) [Erratum- 
ibid. B 223, 544 (1983)]. 

[32] G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, 
M. Sher and J. P. Silva, Phys. Rept. 516 (2012) 1. 




